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Abstract 

As a first step in the aim of arriving at a differential description of neutral Higgs boson 
decays into heavy quarks, h — > QQX, to second order in the QCD coupling as, we have 
computed the hQQ amplitude at the two-loop level in QCD for a general neutral Higgs 
boson which has both scalar and pseudoscalar couplings to quarks. This amplitude is given 
in terms of a scalar and a pseudoscalar vertex form factor, for which we present closed 
analytic expressions in terms of one-dimensional harmonic polylogarithms of maximum 
weight 4. The results hold for arbitrary four-momentum squared, g^, of the Higgs boson 
and of the heavy quark mass, m. Moreover we derive the approximate expressions of these 
form factors near threshold and in the asymptotic regime m^/g^ <^ 1. 
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1 Introduction 



One of the key issues of present-day particle physics is the search for Higgs bosons, that is, 
the search for an answer to the question whether or not the Higgs particle of the standard 
model (SM), or several types of Higgs bosons as predicted by SM extensions, do exist. If Higgs 
particles will be found at the Tevatron and/or at the LHC, then the next important goal will be 
the determination of its, respectively their, couplings to fermions and gauge bosons. This task 
may be achieved to some extent at the LHC, but it is expected that it can be carried out in 
detail only at a high-energy, high- luminosity linear e^e~ collider. There is a vast literature on 
the phenomenology of (non)standard Higgs boson(s) (cf. [1,2] for a recent review), and on the 
interplay of hadron and lepton colliders, in particular as far as Higgs physics is concerned [3] . 

In this paper we are concerned with the couplings of neutral Higgs bosons to quarks (and 
leptons). Apart from the possibility of getting information about the top Yukawa coupling in 
associated production of tih, the Yukawa couplings of neutral Higgs bosons must be determined 
in their decays to 66, cc, and t^t~ (indirect information on the Yukawa couplings can also be 
obtained from the branching fractions of the channels h — > 77, gg). If heavy non-standard 
Higgs bosons exist, experimental investigation of the channel h ^ ti may also be feasible. The 
fermionic decay rates have been calculated quite precisely, both for the SM and for some of its 
extensions, notably the minimal supersymmetric extension (MSSM). As far as QCD corrections 
to the decay rates of scalar and pseudoscalar Higgs bosons into heavy quarks, h QQ, are 
concerned, the order as contributions were determined a long time ago [4-7], and the order 
Og corrections were computed, employing quark mass expansions to various orders, by [8-15]. 
(The hadronic decay rate, i.e., with the gluon channel included, is known to order q;| for ruq — 
{q 7^ t) and ruh <^ mt, see [16,17]). Results on closed quark loop insertions to the scalar [18,19] 
and pseudoscalar [19] vertex functions at two loops are known in the literature already. The 
SM electrowcak corrections to /i — * // were determined in [20-26] and the radiative corrections 
that arise within the MSSM in [27-30]. 

However, not only the partial decay rates but also differential distributions and, in the case 
of tau leptons and top quarks, also the spin polarizations and correlations of the final state 
fermion-antifermion pairs should be predicted as precisely as possible, since these observables 
contain important information about the properties of the decaying boson. For instance, neutral 
Higgs boson(s) may exist that are CP-violating admixtures of scalar and pseudoscalar states. 
In order to find out whether h has quantum numbers J^'^ — 0++, 0~+, or has undefined CP 
parity one may employ suitable spin-spin correlations and polarization asymmetries which were 
proposed and computed, including radiative corrections, in [31-33]. These apply, for light and 
intermediate mass Higgs bosons, to h ^ t~^t~ , while in the case of heavy Higgs bosons the 
channel h ^ tt becomes also relevant. 

We envisage a completely differential description of the (on-shell or off-shell) decays of 
neutral Higgs bosons to heavy quarks, h — > QQX, to second order in the QCD coupling as, 
where Q = c,b,t (this includes the resonant production of Higgs bosons at a future muon 
collider, fi~^fi~ —^h^ QQX). Such a differential description is necessary, as already indicated 
above, for the computation of many observables, especially if phase space cuts are to be taken 
into account. To this order in the perturbative expansion, the amplitudes for h — > QQ, h — > 
QQg, and h QQij, where ij — gg, qq, QQ are required within the SM - and, in addition, 
a suitable jet calculus, which is yet to be worked out for massive quarks at next-to-next-to- 
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Figure 1: Tree-level and one-loop diagrams that contribute to the heavy-quark scalar and 
pseudoscalar form factors The curly line represents a gluon, the double-line denotes the 
heavy quark of mass m. The external quarks are on their mass- shell: Pi = P2 = 



leading order (NNLO) in the QCD coupling (for the case of massless quarks a NNLO subtraction 
method was recently developed [34], cf. also references cited therein). We consider the most 
general neutral Higgs boson state, i.e., a state which can couple both to scalar and pseudoscalar 
fermion currents. This contains scalar and pseudoscalar Higgs particles as special cases. As a 
first step in this project we present in this paper the h — > QQ amplitudes to order a| with 
full quark mass dependence for arbitrary squared four-momentum of the spin zero particle. Of 
course this yields as a special case also the second order QED corrections for the amplitudes 
h —>■ f f, f = quark or lepton. 

The paper is organized as follows: In Section 2 we define the scalar and pseudoscalar vertex 
amplitudes for on-shell Q, Q states. These amplitudes are determined by one form factor each. 
Then we specify the renormalization scheme employed in this paper. In Sections 3 and 4 we 
give the renormalized one- and two-loop scalar and pseudoscalar form factors in the spacelike 
and in the timelike regions, respectively. In Section 5 and 6 we perform the threshold and 
asymptotic expansions, respectively, of the form factors. We conclude in Section 7. 



2 The Higgs-Quark-Antiquark Amplitude 

We consider a general neutral Higgs boson h that couples both to scalar and pseudoscalar 
fermion currents. This includes scalar or pseudoscalar Higgs particles as special cases. The 
Yukawa couplings of /i to a massive quark Q read: 

Cy = [aQQQ + bqQi-f^Q] h, (1) 

where Q and h denote bare fields, mo is the bare mass of Q, V = (v^Gf)-^/^ is the SM 
Higgs vacuum expectation value with Gp being Fermi's constant, and ag, bq are dimensionless 
"reduced" Yukawa couplings whose values depend on the parameters of the specific model under 
consideration. 
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Figure 2: Two-loop vertex diagrams involved in the calculation of the form factors at order 
OicPh). The notation is as in Fig. 1. The single straight and dotted lines represent massless 
quarks and ghosts, respectively. The external quarks are on their mass-shell, Pi = pi = ■ 
Crossed diagrams are not drawn. 
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2.1 Scalar and Pseudoscalar Form Factors 



We consider the decay h{q) — > Q{pi) + Q(j>2) and aim at computing the scalar and pseudoscalar 
parts of the hQQ amplitude, for on-shell Q, Q and for arbitrary squared four-momentum = 
■s = (Pi + ^2)^, to second order in the QCD coupling as and to lowest order in the Yukawa 
couplings. For definiteness we work in QCD with Nf massless flavors and one massive quark 
Q. The interaction implies that the hQQ amplitude V is determined by two form factors: 

KiC2(Pl,P2) = iUciiPl) Tcical?) VC2{P2) , (2) 
TTl / \ 

Tcicalg) = - — (aQFs{s) +ibQFp{s) -f5jSc,c2 ■ (3) 

Here Fs and Fp denote the renormalized scalar and pseudoscalar form factors (whose tree-level 
values are equal to one), m is the on-shell mass of Q, and Ci, C2 are color labels. The expansion 
of these form factors in powers of as/{2n), where as is the standard MS coupling in Nf + 1 
flavor QCD defined at the renormalization scale /i, reads 

We regularize both the ultraviolet (UV) and infrared divergences using dimensional regulariza- 
tion in D = 4 — 2e dimensions. We use a renormalization scheme where the wave function and 
the mass of Q are defined in the on-shell scheme, see below. The UV-renormalized form factors 
dU still depend on e because they are infrared divergent. 

The contributions to F are, for the renormalization scheme employed here, the one-particle- 
irreducible vertex functions shown in Figs. Q and |2l and counterterm diagrams which are not 
displayed. We computed all of them in Feynman gauge. 

The contribution to each form factor can be extracted from V"ciC2 by means of a suitable 
projection. In particular one has, 

Tr [ Id {pi +m)T {p2-m)] , (5) 



2 m (s — Am?) 
iv 



2mbo s 



Tr [ 75 ipi + m) T {p2 - m) ] , (6) 



where p = p^7^ and Id is the unit matrix in Dirac space. 

In what follows Nf denotes the number of light quarks running in the loops Fig. |21 (h), and 
Cf = (N^ - l)/(2iVc), Ca = Nc, Tr = 1/2, where N^ is the number of colors. Notice that the 
scalar and pseudoscalar triangle diagrams analogous to Fig. |21(f), with light quarks circulating 
in the triangles, give zero contribution because a trace over an odd number of 7 matrices is 
involved if the masses of these quarks are neglected in the triangle. 



2.2 Renormalization 

Before specifying our renormalization scheme, a remark on the prescription for the matrix 75 
within dimensional regularization is in order. We use an anticommuting 75 in D dimensions 
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with 7I = 1. This is known to be an adequate prescription in the case of diagrams containing 
only open fermion hnes, because it does not lead to spurious anomalies which must be canceled 
by performing an additional renormalization of the respective current, which in our case is 
the pseudoscalar current. As to the diagrams Figs. [T] and an exception to this rule is the 
pseudoscalar triangle diagram Fig. |21(f). However, this diagram is both UV- and IR-finite and 
may therefore be evaluated in 4 dimensions. It was computed already in [35]. 

We use renormalized perturbation theory with as = Qg/l^ir) being defined as the MS 
coupling in Nf + 1 flavor QCD, while we define the mass m and the wave-function of the heavy 
quark Q in the on-shell (OS) scheme. 

For the renormalization of the hQQ vertex function to two-loop order we need the coupling 
renormalization and the gluon wave function Zf^^ to one- loop, the latter one in the Feyn- 
man gauge. These renormalization constants are well-known and will not be reproduced here 
(c.f., e.g., [36]). 

The renormalization constants of the heavy-quark mass and wave-function are defined in the 
OS-scheme. We need these constants to two-loop order; they were computed in [37, 38] . Using 
the results of [38], we express the one- and two- loop contributions in terms of the renormalized 
MS coupling as- The quark wave function renormalization constant Z^"^ expressed in this way 
may be found in [36]. Here we give only the mass renormalization constant to two-loops 
(mo = Z^'m): We have Z^' = Z^^^s + <05 where 



and 
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(7) 



(9) 



is an overall factor associated with each loop integral in D dimensions. Here and in the following 
the Riemann zeta function of integer argument is denoted by Cn- 
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Further, we need the renormahzation constant Zip for the QQ-gluon vertex to one loop. 

This constant is fixed by the Slavnov- Taylor identity Zip = Z^^ \J~Z^. 

Finally, we come to the overall renormahzation of the hQQ vertex function. Because of the 
interaction the renormahzation of the Yukawa couplings is fixed by the quark mass and 
wave function renormahzation constants. Moreover, as stated above, in our prescription for 75 
no additional counterterms are needed for the pseudoscalar vertex to order a|. Thus, the bare 
and renormalized hQQ vertex functions are related by 

r = z^^z^^Vb. (10) 

The counterterm diagram contributions to F far as their topological structure is 

concerned, identical to those given in [36], where the vector and axial vector parts of the 
Z* — > QQ amplitude to order q;| were computed, and will therefore not be exhibited here. 



3 Renormalized Form Factors for Spacelike s < 

The renormalized form factors at the one- and two-loop level -F/^''' and F^^^'-* [i = S, P) are 
computed with the technique that was applied in [35, 36, 39] to the computation of the two- 
loop photon vertex in QED and to the two-loop vector and axial vector vertices in QCD, 
respectively. Performing the D-dimensional 7 algebra using FORM [40], we obtain the form 
factors expressed in terms of scalar integrals, which in the two-loop case amounts to several 
hundreds. They can be expressed as a combination of only 35 master integrals (Mis) by 
means of the by-now standard Laporta algorithm [41], which makes use of integration-by-parts 
identities [42] , Lorentz invariance [43] and general symmetry relations. This reduction algorithm 
is performed exactly in D = [4 — 2e) dimensions [46] . Once the expression in terms of the Mis 
is found, one can expand the result in powers of e around e = (Z) = 4) and use the values of 
the Mis given in [44,45], which were evaluated with the differential equations method [47]. The 
expressions are given as a Laurent series in e where both UV- and IR divergences, regularized 
by the same parameter D, show up as poles in e. After adding the counterterm diagrams, we 
obtain -F/^'^ and -F/^'^ and the remaining poles in e are due to IR divergences. 

Notice that the sum of the diagrams in Figs. ^ and El and of the counterterm diagrams has a 
non-trivial dependence on the renormahzation scale /i: the counterterms with renormahzation 
constants defined in the OS-scheme carry an overall factor (/i^/m^)^*^; while those with renor- 
mahzation constants defined in the in the MS-scheme carry just one factor {jj? /m?y, coming 
from the loop integration. Due to the mismatch between the e-dependence of the constants in 
the OS-scheme and of those in the MS-scheme, the e-expansion of the result generates terms 
proportional to In (/i^/m^) and In^ (/x^/m^). In this section we give the one- and two-loop 
renormalized form factors choosing the renormahzation scale at the value of the heavy quark 
mass, /i = m. For the terms proportional to In {^j? /m?) and In^ (/i^/m^) we refer the reader to 
Section 13.31 

The one-loop form factors are presented below, including the term of order e while the 
two-loop contributions are given to order (e)°. They are expressed in terms of 1-dimensional 
harmonic polylogarithms H(a; x) up to weight 4 [48,49], which are functions of the dimensionless 
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variable x defined by 

X = = —= . (11) 

V — s + + y— s 

We give our results firstly in the kinematical region in which = s is spacelike (0 < x < 1), 
where the form factors are real. In Section H] we shall perform the analytical continuation to 
the physical region above threshold, s > 4m^, — 1 < x < 0, and explicitly decompose the form 
factors into real and imaginary parts. 



3.1 One-Loop Renormalized Form Factors 

The one-loop renormalized form factors are obtained by adding the contributions of diagram 
of Fig. ^ (b) and its counterterm. 



3.1.1 The Scalar Case 

By setting fi = m and defining 



one finds: 



ao= - z 5 + 



1 — / 1 — 



2(1 + x^) Hj-l, 0, x) (1 + x^) HjO, 0, x) 
^ l-x2 l-x2 ' ^ ^ 



2 (-l + x2-2xC2 + C3 + a;\; 



3) 



1 — X^ 
2 (l + x^) C2H{-l,x) 
1 — x^ 

-1 + 6X-X2 + C2+X2C2) i^(0,x) 



+ 



1 -X2 

8xH{-l,0,x) ^ 4x^r(0,0,x) 



1 — x^ 1 — x^ 

4(1 + x^) i^(-l, -1, 0, x) 2(1 + x^) H{-1, 0, 0, x) 



(12) 



. (l + x^) i^(0,x) 
a-i= -1 ^—^^ ; (13) 

l-x2-C2-x2C2\ 4xi7(0,x) 



1 — x^ 1 — x^ 

2 (l + x^) g(0,-l,0,x) _ (1 +x2) g(0,0,0,x) 

+ 1-X2 1-X2 • 
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3.1.2 The Pseudoscalar Case 

Using an analogous expansion for the pseudoscalar form factor, 



Fi}'\s,e) = C{e) Cf( -a^i + do + e ] , (16) 



one finds: 



a-1 = a-1 ; (17) 
AxH{0,x) 



ao = ao- 



1 — 



-4 (1 + x^) H{-1, -1, 0, x) 2 (1 + x^) H{-1, 0, 0, x) 



ai = ; ^ h 



1 — 1 — x^ 

^2 (1 + x^) H{0, -1, 0, x) (1 + x^) H{0, 0, 0, x) 



1 — x^ 1 — x^ 

2 (l + x2) H{-l,x)C2 



1 — x^ 

i7(0, x) (-1 - 2 X - x2 + C2 + x2 C2) 



1 — x' 



2 



2(-l+X^+C3 + X^ 

+ • ^^^^ 

3.2 Two-Loop Renormalized Form Factors 

The two-loop renormalized scalar and pseudoscalar form factors are obtained by adding the 
contributions of the diagrams in Fig. El (a)-(j), of the crossed ones, and of their counterterms. 
The pseudoscalar contribution of Fig. |21 (f) was already computed in [35] and it is added here 
for completeness. 

3.2.1 The Scalar Case 

By setting fi = m and defining 

i^^'^(6,3) = C2(6){c|(^ + ^+ 6o)+C^C^(^ + ^+ CO 

+CFTnNj[ ^ + ^ + do) + C^T^(^ + ^ + eo 
+0{e) } , (20) 



one finds: 
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1 ^ (l + x^) H{0;x) ^ + H{0,0;x) 



1-x^ 
-1 +^2 + (2 + ^2(2 



{1-x^y 



l-x^ 

(-l + 4x-4x=^ + x^ + C2 + 2^2 (2 + ^^(2) ^(0;^) 

(1 - x^f 

2 (l+x^) H{-l,0;x) {-l + 8x + 8x^ + x'^) H{0,0;x) 



4 (1 + x'^f H{-1, 0, 0; x) 2 (1 + x'^f H(0, -1, 0; x) 



{1-x^y 



(i-x'^y 



3{l+x^f H{0,0,0;x) 



{1-x^y 
1 



20{i-xy{i + xy 



145 + 145 a; - 290 x' - 290 a;^ + 145 x^ 



+145 x^ + 20 C2 - 760 XC2 + 1440 ln(2) x C2 + 520 x'^ (2 
-1440 ln(2) x'^ C2 + 720 x^ C2 - 1440 ln(2) x^ C2 
-540 x^ C2 + 1440 ln(2) x^ C2 + 40 C2 - 118 (2^ 



+ 18 ,T (2'^ - 352 x^ (2^ + 840 x'' C2' + 226 x" (2 
+362 C2^ - 280 Cs - 960 x Cs + 200 x^ Cs 
+200 x^ Cs - 880 x^ C3 - 200 x^ Cs 

(-5 + 19 ,r - 53, r- + (2 //(-!; x) 
(1-x) {l + xf 



3 A 2 



„4 A 2 



3 + 7a;+ lOx^- 10x'^-7x* + 3a;^ 



(l-x)'(l + x)' 

-2 C2 - 10 a; C2 + 47 a;' C2 - 27 x^ C2 - 45 x^ (2 
+5 a;^ C2 + 8 Ca + 24 a; Ca - 56 x^ (s + 56 x^ (s 



-24x^3-8x^(3 



//(0;x) 



4(i + x^)^C3^(i;^) 

(1-X2)' 

(15 -42x + 42x3 - 15x^ + 4(2 + 8x^(2 + 4x^(2) H{-l,0;x) 

(l-x2)' 

2 (-7-7x+16x2-20x^ + 5x^ + 5x5) (2 ^^(0, -1; x) 
(l-x)'(l+x)' 

3 - 12 X + 36 x^ - 22 x^ - 47 x^ + 26 x^ - 3 C2 



1 

(l-x)^(l + x) 
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2 + 2a; + 4 - 4x=^ - + + C2 



+xC2- 12x2(2 + 32x3(2 + 9x^(2 + 13x^C2 H{0,0;x) 
4 

^(l-x)'(l + x)' 
+5 X C2 - 18 x' C2 + 22 x=^ C2 - 3 x^ C2 + a:' C2] ^(1,0; x) 
|4(l + x2) //(-l,-l,0;x) 



1-X2 

(5 + 58x-2x2 + 58x=^ + 9x^) //(-l,0,0;x) 



(1 -X2)" 

2 (3- 22x + 20x2-22x3 + 5x^) H{Q,-l,0;x) 

(1 -X2)2 

(-1 + 31 X + 65 x2 - 9 x3 + 10 x^) ii"(0, 0, 0; x) 

(1 -x)2(l + x)3 
4 (l-6x + 2x2-6x3 + x^) H{0,l,0;x) 
(l-x2)' 

2 (5-16x + 6x2-16x3 + 5x^) H{l,0,0;x) 



(1 - X2)^ 



16(l + x2)2i7(-l,-l,0,0;x) 8(l + x2)'i7(-l,0,-l,0;x) 



(1-X2 



(1 -X2)^ 



12 (1 + x'^f H{-1, 0, 0, 0; x) 4(1 + x^ H{0, -1, -1, 0; x) 



2n2 



(1 -X2)^ 



+ 



(1 -X2)^ 



2 (3 - 5x + 36x2 -16x-^ + 15x^ + 7x5) h{0, -1,0,0; x) 

{l-xf{l + xf 
2 (-l + 3x + 10x2 + 10x=^ + 7x^ + 11x5) i/(0,0,-l,0;x) 

(l-x)'(l+x)=^ 
(7 + 7x + 16x2 + 52x3 + 27x^ + 27x^) //(O, 0, 0, 0; x) 

(l-xfil+x)' 
4 (-1 + 2.V- + 3,ri) //(O. 0. 1. 0; x) 



{1-x^Y 

4 (-1 + ,r) //(O. 1.0, 0: x) 8 (1 + ,7-2)- H{1. 0. -1. 0: x; 



1 + x 



(1-X2)^ 



8 (1 + 3x- 8x2 + 12x3 + h{1,0,0,0; x) 



(l-x)'(l + x)'^ 



8(l+x2)2ij-(l,0, l,0;x) 
(l-x2)2 



11 11 (1 + x2) H{0;x) 

12 ^ 12 (1 -x2) 
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C-1 



49 + 98 - 49 + 18 C2 - 72 x'^ C2 



+54 - 18 Cs - 36 Cs - 18 x^s 



36(1 -x2)2 
2 " 

(1 + x^) (67 - 67 a:^ - 18 C2 + 54 C2) ^(0; a;) 
36(1 

(1 + if(-l, 0; x) 2 ii"(0, 0; x) 



1 — x"^ 1 — x"^ 

(1 + i7(l, 0; x) (1 + x^f H{0, -1, 0; x) 



+ 



l-x^ (1 - x\ 

2 (x2 + x^) H{0, 0, 0; x) (1 + x^f H{0, 1, 0; x) 



{1-x^y 
1 



Co 



540(l-a;)^(l + x)^ 



{l-x^y 

4345 - 4345 x + 8690 x'^ + 8690 x^ - 4345 



-4345 x^ + 2730 C2 + 16860 x (2 - 19440 ln(2) a; C2 
-18900x^2 + 19440 ln(2)a;2C2 - 15840 r'^C2 
+ 19440 ln(2) x^ C2 + 16170 x^ C2 - 19440 ln(2) x^ C2 
-1020 C2 - 81 C2^ - 81 a; (2^ - 3132 a;^ (2^ - 864 a;^ (2^ 
-1917a;^C2^ - 1917x''' (2^ + 4410 Cs - 1530 a; Cs 
+2700 a;^ Cs + 2700 x^ Cs - 7110 x^ Cs - 1170 x^ Cs 
1 

6(1 

+7 C2 + 12 Cs + 24 x^ Cs + 12 Cs 



25 C2 + 162 X C2 - 306 C2 + 162 x^ C2 
H{-l;x) 



1 



121 + 1133 X + 1254 - 1254 x^ - 1133 



54(l-x)^(l + a;)^ 
+ 121 x^ - 99 C2 + 117 X C2 + 999 x'^ (2 + 729 x^ C2 
-1764 x^ C2 + 18 x^ C2 + 351 Cs + 351 x Cs + 594 x^ Cs 



-702 x^ Cs - 405 X* Cs - 405 x^ Cs 



ii'(0;x) 



(1 + x^) (-C2 + x^C2 + 2Cs + 2x^Cs) H{l;x) 
{1-x^f 

(-53 + 102X- 102x3 + 53x^-72 C2- 72x2 C2) ii'(-l,0;x) 

18(l-x2)2 

(-5-5x + 20x2- 16x3 + 7x^ + 7x5) C2^^(0,-l;x) 
(1 -x)^(l + x)^ 

- 67 + 74x + 78x2 -24x3 -11x^-50x5 



18(l-x)^(l+x)' 



+I8C2 + I8XC2 + 54x2C2 + 18x3 C2 + 18x^C2 



(25) 
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+18x^(2 



H{0,0;x) 



{l + x^)\2H{0,l;x) 
{l-x'^f 

2 jl + x^) -4(2 + 2x^(2) H{l,0;x) 

{1-x^f 

52 + H{-l,-l,0;x) 



3 (1 - ,t2) 

(61 + 175 X - 35 x2 + 79 a;^) i/(-l, 0, 0; x) 

6 (l-x) {1 + xf 
6 (l+a;2) H{-l,l,0;x) 
1-^2 

2 (1 + 10X + 25^2 + 34x3) iJ(0,-l,0;a;) 

3 (l-x) (1 +x)^ 
(11 - 2 X - 53 x^ - 74 x^ + 62 x^) H{0, 0, 0; x) 
6 (1 -x) (1 + x)^ 

4x {1 + .V + 2.V-) //((). 1.0:, 7:) 



(l-x) (1 + x)' 
6 (l + x^) //(l,-l,0;x) 4 (-3 + x^) H{1, 0, 0; x) 



1 — x^ 1 — x^ 

2(1 + x^) i7(l, 1, 0; x) 4(1 + x^)^ i7(-l, 0, -1, 0; x) 



1-X2 



(1 -X2 



2 (3 + 4x2 + x^) H{-1, 0, 0, 0; x) 4 (1 + x^ H{-1, 0, 1, 0; x) 



(1 -X2)^ 



(1-X2 



10(l + x2)^iJ(0,-l,-l,0;x) 

(1 - ,t2)2 

2 (1 + x + x2 + 15 ;r'^ + 7x4 + 7x5) i/(0, -1, 0, 0; x) 



(l-x)^(l + x)^ 



6(l + x2)^ii"(0,-l,l,0;x) 
(l-x2)' 

2 (1 + X + 22 x2 + 2 x^ + 11 x^ + 11 x^) H{Q, 0, -1, 0; x) 

(1 -x)^(l + x)^ 
3x2 (1 + 7a; + 4x2 + 4x3) h{Q,Q,Q,Q; x) 



(1 -x)^(l + x)' 



2 (1 + 8x2 + 7x4) i7(0,0,l,0;x) 6 (1 + x2)^ /7(0, 1, -1, 0; x) 



(1 -x2)^ 



+ 



(1 -X2)^ 



2 (3 + 3 X - 8 x2 + 16 x^ + x^ + x^) if(0, 1, 0, 0; x 
(l-x)^(l + x)^ 
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2 (1 + x^Y H{0, 1, 1, 0; x) 4 (1 + x^ f H{1, 0, -1, 0; x 



+ 



{1-xY {1-xY 
2 (-5 -4x2 + i/(l,0,0,0;x) 4 (1 + x^)^ 0, 1, 0; x) 



(1 -x' 



2^2 



(1 -X 



2^2 



_1 _ (1 + x2) ij'(0;x) 
~3 3 (1 - x2) ' 
5 5 (l + x^) i/(0;x) 
9^ 9 (1 - x2) ' 
-49 + 49 x^ - 6 C2 - 72 X C2 + 66 x^ (2 + 36 (3 + 36 x^ (3 



27 (1 - x2) 

4 (l + x^) C2ij'(-l;a:) 

3(1- x2) 

2 (14 -96x+ 14x^ + 9(2 + 92:2(2) i^(0;x) 
27 (1 -x2) 

4 (5- 12x + 5x2) i/(-l,0;x) 2 (5- 12x + 5x2) H{0,0;x) 



+ 



9(l-x2) 9(l-x2) 
(l + x^) H{-l,-l,Q;x) 4 (l + x2) i/(-l,0,0;x) 



3 (1 - x2) 3 (1 - x2) 

4(1 + ,/•-') HiO. -1. 0: ,r) 2(1 + .1'-) H(0. 0. 0: ,t) 



3 (1 - x2) 



3 (1 - x2) 



0; 
0; 

J [ - 2035 - 1435 X + 5505 x^ + 4905 x^ 

135(1 -xf (1 +x) L 

-4905 x^ - 5505 x^ + 1435 x^ + 2035 x^ + 810 (2 + 270 x (2 
-9810 x^ C2 + 9450 x^ C2 - 810 x^ (2 + 5490 x^ C2 
-4590 x^ C2 - 810 x^ C2 + 432 x (2^ + 2160 x^ 
+2592 x^ C2^ + 2592 x^ (2^ + 2160 x^ (2^ + 432 x^ (2^ 
+8640 x^ Ca + 8640 x^ Ca - 8640 x^ Cs - 8640 x^ Ca 
2 

^27(1 -x2)=^ 

+9(2-9x2(2-9x^(2 + 9x^(2 - 162 X (3 - 648x2 (3 



28 + 144 X + 212 x2 - 768 x^ + 212 x^ + 144 x^ + 28 x 



+756 x^ (3 - 648 x^ (3 - 162 x^ (3 
32xiJ(-l,0;x) 



H{Q;x) 



l-x2 
2 



9(l-x)'(l + x)' 



5 - 23 X + 101 x2 + 531 x-^ - 27x^ 
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-353 - 229 x^ + 5x'^ + 18 x (2 + 90 x^ (2 + 468 x^ (2 



+468x^(2 + 90x^(2 + 18x^(2 
8xH{0,l;x) 



H{0,0;x) 



1 -X2 



8x (2 + 4x + 2 ,x2 + (2 + 6 ,T (2 + x2 (2) H{1, 0; x) 

(1 - x) (l + x f 
128 x^ H{0, -1, 0; x) 2 (-1 + 48 x^ + x*) H{Q, 0, 0; x) 



(l-x2)^ 



3(1 -x2 



64x2if(0,l,0;x) 128,x2i7(l,0,0;x) 256x3 iJ(0, 0, -1, 0; x) 



(1 -X2)^ 



+ 



(1-X2)^ 



(1-X2) = 



2x (1 + x) if(0,0,0,0;x) 12x (1+x) H{0,0,0,l;x) 



(1-xr 

8x (l + 6x + x2) i7(0,0,l,0;x) 



(1-xr 



(l-x) (l + x)-^ 



8x (l + 4x-26x2 + 4x=^ + x^) //(0, 1, 0, 0; x) 



(l-x2f 
8x (l + 6x + x2) //(l,0,0,0;x) 
(l-x)(l + x)^ 

3.2.2 The Pseudoscalar Case 

As done in the scalar case, by choosing /i — m and defining 



(32) 



F^;\e,s) = C\e){cU^-^ + — + h] + CfCa( ^ + — + co 



C-2 , C_i 



+C^r^7V,('^ + ^ + do) + C^T^r^ + ^ + eo 



+0(e) , 



one finds: 



(33) 



6_2 = &-2 ; 



4xi/(0;x) 8x (1 + x2) i/(0,0;x) 
4 (l + x2) iJ(-l,-l,0;x) 



1 — x^ 

(5 + 14x + 46x2 + 14x3 + 9x^) iJ(-l,0,0;x) 
(l-x2)2 



(34) 
(35) 
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2 {3-Ux-4x^ - Ux^ + 5x^) H{0,-l,0]x) 

{1-x^f 

(1 - 9 a; - 9 + 71 - 16 + 10 x^) H{0, 0, 0; x) 
{l-xf{l + xf 

4 (l + a;2) H{0,l,0;x) 2 (5 - 12a; - ISx^ - 12x=* + 5x^) H{l,0,0;x) 



(l + xY 



+ 



(1-^2)^ 



_^16{l + x^f H {-1,-1, 0,0; x) ^8{l + xy H{-l,Q,-l,Q;x) 



(1-^2)^ 



{1-x^y 



12 (1 + x^f H{-1, 0, 0, 0; x) 4(1 + x^ H{0, -1, -1, 0; x 



2^2 



(l-a;2)^ 



+ 



(1-^2)^ 



2 (-3 + llx-12x2 + 8x3 + x^ + 7a;5) i/(0, -1, 0, 0; x) 

(1 - x)^ (1 + a;)^ 
2 (l-Sx-lOx^ + lOx^-lSx^ + llx^) //(0,0,-l,0;a;) 

(1-xf (1 + x)' 
(-7 + 7a: - 24^2 + 44a;=^ - 27x^ + 27 a;^) H{0, 0, 0, 0; x) 

{l-xf{l + xf 
4 (-1 + 2 a;2 + 3 x^) H{0, 0,1,0; x) 



{1-x^f 

4 (l-4a; + 2,T2-4x'^ + x^) H{0,l,0,0;x) 



{l-xy {1 + x) 
8{l+x^f H{l,0,-l,0;x) 
{1-x^f 

8 (1 + 2 X + 2 x2 - 2 + x^) H{1, 0, 0, 0; x) 

(1 -a:2)2 

8(l + a;2)^ii'(l,0,l,0;x) 

(1 - x^f 

(5- 18a; + 54;7;2 - ISx-'^ + a;^) H{-l;x)C2 

(1 - X2)2 

2 (7 - 7a; + 8a;2 + 4,r'^- 5x^ + 5x5) i/(0, -1; a;) C2 

(l-x)^(l + a;)^ 
4g(l,0;a;) (-2 + 2a;^ + C2 + 4a: C2 + 2a;2 (2 - 4a;3 C2 + a;^ C2) 

(l-a;2)2 

//(-I, 0; a;) (15 + 2 a; - 2 a;=* - 15 a;^ + 4 C2 + 8 a;2 C2 + 4 x^ C2) 

(l-a;2)2 

- 3 + 14a; + 6a;2 - 8a;3 - 35a;^ + 26a;^ + 3C2 



(l-a;f (l + a;f 



-7a;C2 + 4a;2 C2 + 24 a;^ (2 - 17 a;^ C2 + 13 a;^ C2 



H{0,0;x) 
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A{l + x^fH{l;x)C^ 

\ ^ f - 145 + 145 X + 290 x^ - 290 x^ - 145 + 145 x^ 

20(l-a;)^(l + x)^ L 

-20 C2 + 240 xC,2 - 1600 x^ C2 + 1640 x^ C2 - 300 x'^ C2 + 40 x^ C2 
-480 a; ln(2) (2 + 1440 ln(2) C2 - 1440 ln(2) C2 
+480 x^ ln(2) C2 + 118 (2^ - 254 x (2^ + 48 (2^ + 536 x^ (2^ 
-498 x^ C2^ + 362 x^ (2^ + 280 Cs + 280 Cs - 280 x^ Cs 
-80 x^ Cs - 200 Cs 
H{0;x) 



(l-a;)^(l + a;)^ 



3-x-2x^-2x''-x* + 3x^ + 2C2-6xC2 



-11 C2 + 73a;3 C2 - 15a:X2 + 5a;^ C2 - 8C3 - 8a;C3 



(36) 



c-2 = C-2 ; 



C-1 = c_i ; 



(37) 
(38) 



Co 



-52 (l + a;2) i7(-l, -1, 0; a;) 
3 (l-a;2) 

(-61 - 6 X + 66 ^2 - 6 a;3 + 79 x^) (-1, 0, 0; a;) 

6(l-a;2)' 
6 (l + a;2) //(-1, 1,0; a;) 

_2 (-1 - 15a;- 15a;2 - 15a;3 + 34a;4) H{0,-l,0;x) 

3(l-a;2)' 

(11 - 35 a; + 45 a;2 + 189 a;^ - 128 x^ + 62 a;^) H{0, 0, 0; a;) 

6 (1 - x)^ (1 + a;)^ 
4x (l + a; + 2a;2) H{0,l,Q;x) 



(i-x) {i + xy 

6 (l + x2) iJ(l,-l,0;a;) 

4 (3-2a--4a-2-2a;3 + a;^) if(l,0,0;a;) 

(l-a;2)' 
2 (l + a;2) i/(l,l,0;x) 



l-a;2 



A {1 + x'^f H{-1, 0,-1,0; x) 
(l-a;2)' 
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2 (3 + 4a;2 + x^) H{-l,0,0,0]x) 
4(1 + a;2)^i/(-l, 0,1,0; 
10(l + a;2)^i7(0,-l,-l,0;x) 

2 {-l + x-5x^ + llx^ -7x^ + 7 x'') H{0, -1,0,0; x) 

{l-xf{l + xf 
6{l+x^f H {0,-1, 1,0; x) 
{1-x^f 

2 {-1 + x - Ux'^ + lOx^ - llx"^ + llx^) H{0,0,-l,0;x) 

{l-xf{l+xf 
x^ (-7 + 17x- 12x2 + 12x3) H{0,0,0,0;x) 

{l-x)^{l + xf 
2 {l + Sx^ + 7x^) g(0,0,l,0;a:) 

{1-x^f 
6{l+x'^f H(0, 1,-1,0; x) 
{1-x^f 

2 {-3 + 3x + 8x^ -x'^ + x^) H{0,l,0,0;x) 

{l-xfil + xf 
2(1 + A-'^f H(0. 1. 1.0: x) 

{1-x^f 
4(1 + .v'^ f H{1. 0. -1.0; x) 

{1-x^f 
2 (-5 - 4 x^ + x^) H{1, 0, 0, 0; x) 

(l-x'^f 
4(l+x2)^/j'(l,0,l,0;x) 
(l-x^f 

{5-5x + 4:X^ + 8x^ -7x^ + 7 x^) H{0, -l;x) C2 
(l-a;)'(l + a;)' 

(l + x^f i/(0,l;x)C2 

{1-x^f 

H{-l,0;x) (-53 - 54 x + 54 a;3 + 53x^- 72 (2 - 72x2(2) 

18(l-x2)2 
2 (l + x2) ij"(l,0;x) (l-x2 -4(2 + 2x^(2) 
(l-x2)2 
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H{0,0;x) 



18(1-,t)''(1 + x)^ 



67 - 148 X + 36 -90x^ + 185 - 50 - 



+18 xC2- 126 C2 - 54 C2 - 18 x^ (2 + 18 (2 
{l + x')H{l;x){-C2 + x'C2 + 2C3 + 2x'Cs) 



H{-l-x) 



{l-x-'f 

25C2 + 54xC2- 162x'C2 + 54x^C2 + 7x^C2 



6 {1-x^Y 
+12(3 + 24x2(3 + 12x^(3 

\ . [4345 - 4345 X - 8690 x^ + 8690 x"^ + 4345 x^ 

hm{l-xf{l + xf L 

-4345 - 2730 (2 - 1320 a; (2 + 25200 x^ (2 - 23580 x^ (2 

+3450 (2 - 1020 x^ (2 + 6480 x ln(2) (2 - 19440 x^ ln(2) (2 

+19440 x^ ln(2) (2 - 6480 x^ ln(2) (2 + 81 (2^ - 81 a; (2^ 

-108x2(2^-4104x3(22 + 1917x^(2^ - 1917x^(2^ - 4410(3 

+4950 X (3 - 2700 x^ (3 + 2700 x^ (3 + 630 x^ (3 - 1170 x^ (3 

//(0;x) 
54(1 -x)^(l + x) 
+99 (2 - 315 X (2 - 135 x2 (2 + 2241 x^ (2 - 612 x^ (2 
+18 x^ (2 - 351 (3 + 351 X (3 - 162 x^ (3 - 270 x=^ (3 

+405x^(3-405x^(3] ; 



121 - 283 X + 162 x' + 162 x^ - 283 x^ + 121 



d. 



Jo 



= d-i ; 



-20 (l + x2) //(-l,0;x) 10 (l + x2) /j-(0,0;x) 
9 (l-x2) ^ 9 (l-x2) 

8 (1 + x2) H{-1, -1, 0; x) 4 (1 + x^) H{-1, 0, 0; x) 



3 (1 -x2) 



3 (l-x2) 



4(1 + x2) iJ(0, -1, 0; x) 2(1 + x2) i/(0, 0, 0; x) 



3 (1 -x2) 



+ 



3 (l-x2) 



4 (l + x2) ij'(-l;x)(2 2 (1 + x2) H{Q- x) (14 + 9 (2) 



3(l-x2) 27(l-x2) 
-49 + 49 x2 - 6 (2 + 66 x2 (2 + 36 (3 + 36 x2 (3 
27 (1 - x2) 



= e_2 = 0; 
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e_i = e_i = ; 



(44) 



eo 



2 {l-2x + 2Qx^ -2x^ + x^) H{0,0,0;x) 2x H{0,0,0,0;x) 



3(l-x)'^ (1 + a;) 



1 -x2 



12a;//(0,0,0, l;x) 8 x i7(0, 0, 1, 0; x) 8a;i/(0, 1, 0, 0; x) 



1 -a;2 



1 -x2 



1 -x2 



8xi/(l,0,0,0;x) 8 2;/7(l,0;x)C2 



1 — 

2if(0,0;x) 



1 — x^ 

5 + 22 X + 59 + 116 x^ + 59 x^ + 22 x^ + 5 x^ 



9(l-x)'(l + x)' 

-18 X C2 - 36 x^ C2 + 36 x^ C2 + 18 x^ C2 

2/7(0;x) 
27(l-x2)3 

+207 x' C2 + 432 x^ C2 + 207 x^ C2 + 9 x^ C2 - 162 x (3 
+324x^(3 - 162x^C3 



28 + 24 X - 76 x^ + 48 x^ - 76 x^ + 24 x^ + 28 x^ + 9 C2 



135(1 -x)^(l + x) 



2035 + 4190 X - 2035 x^ - 8380 x^ - 2035 x^ 



+4190 x^ + 2035 x^ - 810 C2 - 1080 x C2 + 4410 x^ C2 + 12240 x^ C2 
+4410 x^ C2 - 1080 x^ C2 - 810 x^ C2 - 432 x (2^ - 864 x'^ (2^ 

+864x^2^ + 432x^(2^^ 



3.3 Form Factors for fi ^ m 



(45) 



At the one-loop level, the dependence of the form factors on the renormalization scale /z appears 
in an overall factor (/i^/m^)^: 



C{e) 



^,P) 



(46) 



with the expression of Fg^''\e,s) and Fp^\e,s) given in Eq. (fT^ and Eq. (fT^ respectively. 

At the two-loop level, the coupling constant renormalization generates logarithms of the 
ratio (/i^/m^). Factoring an overall (//^/m^)^^, one has [36]: 




-.^F^'%-,s)W(^]\ (> = S.P). (48) 



where the functions F^'^^\e, s) are given in Eq. (j20|l and Eq. (j33|l . and /So is the first coefficient 
of the QCD /5-function, 

^^^ nC.-4T,W^l) ^ (49) 

D 

4 Analytic Continuation above Threshold 

For apphcation to Higgs decay, the form factors must be analytically continued into the timelike 
region, above the QQ threshold, s > 4m^, where they develop imaginary parts. That can be 
made by the substitution s ^ s + ie, which implies 

X — >• —y + ie , (50) 

where 

^ _ y/j - y/g - 4m2 ^^^^ 
^/s + \/ s — Am? 

with y = 1 at s = Am? and ?/ — as s — * oo. 

The real and imaginary parts of the form factors are defined through the relations: 

F,{e,s + it) = ^F,{e,s)+i7rQFi{e,s) {i = S,P). (52) 

In the following two sections we will give the real and imaginary parts of the one- and 
two-loop analytically continued form factors for /i = m. The renormalization scale dependence 
follows from the pattern outlined in Section IH.HI 

4.1 One-Loop Form Factors above Threshold 

In this section we give the analytic continuation of the coefficients a^, corresponding to the 
scalar form factor defined in Eq. (|T^ . and at, corresponding to the pseudoscalar one defined in 
Eq. (HEl). 



4.1.1 The Scalar Case 



^a^,= (53) 



1 + HiO;y) 
l-y^-4(2-^y'(2\ ^yH{0-y) 



1 — ?/2 / ^ — 

l+y^) H{0,0;y) 2 {1 + y^) H{l,0-y) 



X _ y2 1 — 1/^ 

3) 



(54) 



2 (-l+y^ + 8y(2 + (3 + y\: 



i-y' 
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8_(l + 2/2) C2H{l-y) 4yH{0,0-y) 8yH{l,0-y) 



X _ ^2 i — y^ ^ — 

(1 + /) HiO.O.O-y) 2 (l + y2) //((). 1.0;y) 



1 _ y2 1 — 

2 (1 + y2) 0, 0; y) 4(1 + y') H{1, 1, 0; y) 



1 + 2^ 

1-1/2 

^ _ -Ay (1 + y^) /f(0;y) 2(l + y^) H{l-y) 

TT^^ rrp ' ^^^^ 

l + 6y + y^-2C2-2y"C2 
l-y2 

4y//(0;y) 9>yH{l-y) (1 + y^) //(o, 0; y) 



1 — 1/2 1 — 1/^ 1 — 

2 (l + y2) /J(0,l;y) 2 (1 + y^) iJ(l,0;y) 



1 — ?/2 1 — y^ 

4(l+y2) i/(l,l;y) 



l-y2 

4.1.2 The Pseudoscalair Case 



3?ao 



1 _ y2 1 — y^ 

l-y'-4C2-4y2C2 



l-y2 

(1 + y2) i7(0, 0, 0; y) 2(1 + y^) i/(0, 1, 0; y) 



1 _ y2 1 — y^ 

2 (1 + y2) H{1, 0, 0; y) 4(1 + y^) //(1, 1, 0; y) 



1 — y2 1 — y^ 

8(l + y2) //(l;y)C2 



1 -y2 

i/(0;y) (-l + 2y-y^ + 4C2 + 4y^C2) 
l-y2 

2 (-l+y2+^3+y2^3) 

l-y2 



(55) 



$Ja_i = ; (56) 



(58) 



3?a_i= 3?a_i ; (59) 
(l + y2) i/(0,0;y) 2 (l+y^) if(l,0;y) 



(60) 



^t^^^JL^L^I^^^ . (61) 
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53 a_i = 53 a_i ; (62) 



^«o= — -2 T~;r2 ' (^^) 



$5 ai 



(1 + y^) H{Q-y) _ 2 {I + y^) H jl; y) 

1 - ?/2 1 _ y2 

(1 + 2/2) iJ(0,0;y) 2 (l+y2) /7(0,l;y) 



2 (l + y2) i7(l,0;y) 4 (1+2/2) H{l,l-y) 



I — y^ 1 — y"^ 

l-2y + y'-2C2-2yH2 
1-2/2 



(64) 



Our results Eqs. (53), (54) and Eqs. (59), (60) in the scalar and pseudoscalar cases, agree 
with those of [4] and [7], respectively. 

4.2 Two-Loop Form Factors above Threshold 

In this section we give the analytic continuation of the coefficients 6,, q, di, Cj, of the scalar 
form factor defined in Eq. (pUj) and bi, Ci, di, Cj of the pseudoscalar one, defined in Eq. (jSSI)- 

4.2.1 The Scalar Case 

1-22/2 + 2/4-6(2-12 2/2(2-62/^2 



2(1-2/ 

(1 + y') H{0;y) {l + y')'H{0 , 0; y) 
l-y' (1-2/2) 

1 _ 22/2 + 2/4 - 4C2 - 242/ C2 - 24 2/3(2 + 42/^2 



^ — — — + — 7, — 7^ — ; (65) 



(_1_4^ + 4^3^^4^10^^^ 20 2/2C2 + 10|/'C2) H{0-y) 
12 (1 + 2/2)' C2H{l;y) (-1 - 8 2/ - 8 y^ + y^) H{0, 0; y) 



(1-2/2)' (1-2/2)' 
2 (1 + 2/2) H{1, 0; 2/) ^ 3 (1 + 2/')' i/(0, 0, 0; y) 



1-1/' ■ (1-y')' 

2(l+2/2)'g(0,l,0;2/) ^ 4(l+2/2)'g(l,0,0;2/) 

2 + ^T2 ; (66) 



^bn 



(i-yr (1-y') 
1 

~20(l-2/)'(l+2/) 
+ 145 2/^ + 160 C2 - 40 2/ C2 + 1440 ln(2) 2/ C2 + 1640 2/' C; 



- 145 + 145 2/ + 290 2/' - 290 2/^ - 145 y^ 
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+1440 ln(2) C2 + 2040 (2 - 1440 ln(2) C2 
-2280 C2 - 1440 ln(2) y^ C2 - 1520 y^ C2 - 272 (2^ 



+168 1/(2' - 848 y^ C2' + 480 y' C2" - 496 C2 
+392 C2' + 280 C3 - 960 7/ Ca - 200 y^ Ca + 200 Ca 

+880/(3 - 2001/5(3 
2 

-2Ca-4y2Ca-2y^Ca] H{-l-y) 

-J. ^77— ^f-3-7y+10y2 + i0|/3-7|/^-32/5 

+103 y C2 - 148 y2 C2 - 45 y" (2 + 25 + g (3 
-24 y Ca - 56 y^ Ca - 56 y' (s - 24 y^g + 8 y^ C3] H{0; y) 
2 (5-114y-54y2-114y3 + 13y^) C2i^(l;y) 

-2-2y + 4y2 + 4y3-2y^-2y5 



3 A 2 



.4 A 2 



15C2 + 48yC2 + 18y'C2 + 48y=^C2 + 15y'C2 



(i-y)'(i + y)' 



-5C2 + 13yC2 + 30y2C2 + 50y=^C2 + 3y^C2 



+5y'C2 



^(-1,0; y) 



12 (1 + y) C2i^(0,-i;y) 



3 + 12y + 36y2 + 22y3_47y^-26y5 

(1-?/) (1 + ?/) 
-24 C2 + 20 y C2 - 60 y^ C2 + 124 y^ C2 - 72 y^ C2 

+68y5C2] i/(0,0;y) 

8 (-4-2y-23y2-7y=^-10y" + 4y5) C^HjQ.l-y) 

(1-yf (1 + yf 

15 + 42 y - 42 y2 - 15 y^ + 40 C2 + 80 y^ C2 



(1 - V'f 
+40y4C2 



^(1,0; y) 



48(l + y2)^C2i^(l,l;i/) 

2 (5 + 16y + 6y2 + 16y=^ + 5y^) ij-(-l, 0, 0; y) 

(1 - y'f 

4 (l + 6y + 2y2 + 6y=^ + y^) //(0,-l,0;y) 

(i-y^f 
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{l-yf{l+yf 
2 (3 + 22 1/ + 20 + 22 + 5 y^) h{0, 1, 0; y) 

(5 - 58 1/ - 2 1/2 - 58 1/3 + 9 y') H{1, 0, 0; y) 

4{l + y') H{l,l,0;y) 8 {1 + y'f H{^1,0, -1,0- y) 

8 {-l + 3y + 8y^ + 12y^ + y^ + y^) 0, 0, 0; 

(1-l/f (1 + l/f 
8 (1 + y^f H{-1, 0, 1, 0; y) 4(1 + |/) i7(0, -1, 0, 0; y) 



(l-y2)2 l-y 

4(-l + 2y2 + 3y^) ff(0, 0, -1, 0; y) 

(1 - y'f 

(-7 + 7 y - 16 y2 + 52 y3 _ 27 + 27 y^) g(0, 0, 0, 0; y) 

(i-yf(i + yf 

2 (l + 3y-10y2 + 10y3-7y^ + lly^) /J(0, 0, 1, 0; y) 

(l-2/f(l + ?/)' 
2 (-3-5y-36y2-16y3-15y^ + 7y^) i/(0, 1, 0, 0; y) 

{i-yf{i + yf 

4 (1 + y2)2 i/(0, 1, 1, 0; y) ^ 12 (1 + y^)^ H{1, 0, 0, 0; y) 



(1 - y2)2 (1 - y2 

8(l+y2)2g(l,0,l,0;y) 16(l + y2fff(l,l,0,0;y) 

(1 - y'f (1 - y^)' 

11 11 (l + y2) g(0;y) 

12 ^ 12 (1 - y2) 

^ -49 + 98y"-49y^ + 18C2 + 144y^C2 



36(l-y2)2 



-162 - 18 Cs - 36 _ 18 y\3j 

(l + y2) (-67 + 67y2 + i8C2 + 162y2C2) i/(0;y) 
36(l-y2)2 

(l + y2) iJ(-l,0;y) 2y' H(0,0;y) (1 + y^) i/(l,0;y) 



1 — y^ 1 _ y2 1 — y2 

^ (1_+ y2)2 i/(o, -1, 0; y) 2 (y2 + y^) ^(q, q, 0; y) 



(l + y2)'i/(0,l,0;y) 



(1-1/2)' 
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4345 - 4345 y - 8690 + 8690 y^ 



5i0{l-yf{l + yf 
+4345 - 4345 y^ - 8760 C2 + 10200 y C2 - 19440 ln(2) y C2 
+25920 C2 - 19440 ln(2) y^ C2 - 13680 y^ C2 
+19440 ln(2) y^ (2 - 17160 y^ (2 + 19440 ln(2) / C2 
+3480 y^ C2 + 1701 C2' - 1701 y C2' + 10422 C2' 
-19494 y^ C2' + 13257 y^ C2' - 13257 C2' - 4410 Cs 
-1530 y C3 - 2700 y' Cs + 2700 y^ + 7110 y^ Cs 
-1170 y^Ca 

37C2-48y2C2 + lly'C2-2C3-4y'C3 



-2y^C3 



H{-l;y) 
1 



54(l-y)^(l + y)^ 



121 - 1133 y + 1254 y^ + 1254 y^ 



- 1 133 - 121 yS + 198 C2 + 234 y C2 - 378 y^ (2 

-162 y3 C2 + 1908 y^ C2 + 1656 C2 

-351 Cs - 351 y C3 + 594 y' (3 + 702 y^ (3 - 405 y" (3 



+405y^C3 



3(l-y2) 



H{0;y) 

-52C2 + 45yC2 + 81y2C2 + 45y^C2 



+61y^C2 + 3C3 + 6y=^C3 + 3yX3 



2 (1 + y^) (-1 + y2 - 11 + H(-l, 0; y) 

(-19 + 19y + 46y' + 98y3 - 7y^ + 7y5) (2 ij"(0, -1; y) 

{i-yf{i + yf 

--— \ [67 + 74y - 78y^ - 24y=^ + 11 y' - 50y^ 

18(l-y) (1 + y) L 

-18C2 + 18yC2-216y^C2 + 1152y3C2 - 666y^C2 



+666^/(2 



F(aO;y) 



-11 + 11 y + 14 y^ + 106 - 35 y^ + 35 y^) (2 ij"(0, 1; y) 
(l-yf(l + yf 
- 53 - 102 y + 102 y^ + 53 y^ + 252 (2 + 360 y^ C2 



18(l-y2)^ 



+108y^C2 



^^(l,0;y) 
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2 (1 + y^) H{-1, -1, 0; y) 4 (-3 + y^) H{-1, 0, 0; y) 

X _ y2 1 — 

6 jl + y^) if(-l,l,0;j/) Ay {l-y + 2y^) H{0,-l,0;y) 

{l-yf{l+y) 

(11 + 2 y - 53 + 74 + 62 y^) H{Q, 0, 0; y) 

6{l-yf (l + y) 
2 (-l + 10y-25y^ + 34y3) ff(0,l,0;y) 

3(l-yf (l + y) 

6(l + y^) g(l,-l,0;y) 

l-y2 

(-61 + 175 y + 35 y^ + 79 y^) ij-(i, p; y) 

6{l-yf (l + y) 
52 (1 + y^) i/(l, 1, 0; y) 4 (1 + y^)' H{-1, 0, -1, 0; y) 

3(1-1/^) (l-y2)^ 

2(-5-4y^ + y^)g(-l,0,0,0;y) ^ 4 (1 + y^)^ g(-l, 0, 1, 0; y) 

2(l+y^)^g(0,-l,-l,0;y) 
(1-1/^)' 

2 (-3 + 3y + 8y^ + 16y3-y^ + y^) g(0, -1, 0, 0; y) 

(l-y)^(l + y)^ 

6(l + y^)^if(0,-l,l,0;y) 2 (1 + 8y^ + 7y^) g(0, 0, -1, 0; y) 

(i-y'f {1-y'f 

3y' (-l + 7y-4y^ + 4y3) g(P,P,P,P;y) 
(l-y)^(l + y)^ 

2 (-l+y-22y^ + 2y3-lly4 + lly5) g(P,P,l,P;y) 

ii-yf{i + yf 

6{l+y^fH{0,l,-l,Q;y) 

a-y'? 

2 (-l+y-y2 + i5y3-7y4 + 7y5) g(0,l,0,0;y) 

(l-y)^(l + y)^ 
10(l + y^)'iJ(0,l,l,0;y) 4 (1 + y=^)=^ (1, 0, -1, 0; y) 

{i-y'f a^y'f 

2 (3 + 4y^ + y^) if (1, 0, 0, 0; y) 4 (1 + y^)' F(l, 0, 1, 0; y) _ 
{1-y^f [l-y^f 

1 (l + y^)g(0;y) 
3 3(l-y2) ' 
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5 5{l + y^) H{0-y) 

9 9{l~y') 



49 + 49 y2 + 84 C2 + 288 y (2 + 156 (2 



27 

36(3 + 361/2(3 

4 (-7-48|/-7j/2 + 9C2 + 9|/2C2) H{0;y) 
27 (1-1/2) 

16 (1 + 1/2) C2i^(l;i/) , 2{5 + 12y + 5y^)H{0,0;y) 



+ 



3(l-y2) 9(l-y2) 
4(5 + 12y + 5y2)/f(i,0;y) , 2 (1 + y^) //(o, 0, 0; y) 



+ 



9(l-y2) 3(l-y2) 
i{l + y') H{0,l,0;y) , 4 (1 + y^) //(i, Q, 0; y) 



3 (1 - y') 
8 (1 + 1/2) i7(l,l,0;y) 



+ 



3 (1 - y') 



3(l-y2) 



3f?e_2 = ; 
3?e_i = ; 



135 (1 - y^y 



- 2035 - 600 y + 4905 y^ - 4905 y'^ + 600 y^ 



+2035 y^ + 1260 C2 - 1080 y C2 - 19980 y^ (2 + 18360 y^ C2 
+19980 y^ C2 - 17280 y^ C2 - 1260 y^ C2 + 783 y C2' 
-3132 y2 (2^ + 32346 y^ C2' - 3132 y^ C2' + 783 y^ C2' 

+8640 y^Ca- 8640 y^Ca] 
, 384y2C2i/(-l;y) 



(i-y2) 

4 



14 + 72y- 106y2- 384y2- 106y^ + 72y^ 



27(l-y2)^ 

- 14 y6 + 9 C2 - 657 y2 C2 + 639 y^ C2 + 9 y^ C2 - 81 y Ca 



+324 y2 Ca + 378 y=^ Ca + 324 y^ Ca - 81 y' Ca 



^(0;y) 



16y (-l + 2y-y2 + (2-6yC2 + |/2C2) g(-l,0;y) 
(l-y)'(l+y) 

l + 2y2-y4 + 3C2-12yC2-78y2C2 



-12y=^C2 + 3y^C2 i^(0,-l;y) 
2 



5-23y- 101y^ + 531y^ + 27y^- 353y' 



9(l-y)^(l+y)^ 
+229 y6 + 5 y^ - 9 y C2 + 45 y2 C2 + 198 y' C2 - 198 y' C2 



(72) 



(73) 

(74) 
(75) 
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-45^/5(2 + 9/ C2ji^(0,0;y) 
32 y H{1, 0; y) 128 y^ H{-1, 0, 0^ ^ 64 y^ H{0, -1, 0; y) 



(1-1/2 



2 (-l + 48y2 + y4) H{OMy) 128y^ H {0,1,0; y) 



8y (1-6I/ + I/2) Hi-l,0,0,0;y) 

(l-yfil+y) 
8y {l-4y-26y^-4y^ + y^) H{0, -l,0,0;y) 

{1-y'? 

8y {l-6y + y^) H{0,0,-l,0;y) , 12 {-1 + y) y H{0,0,0, -l;y) 



{l-yY {1 + y) 



+ 



i^ + yy 



2 (-1 + y)y H{0, 0, 0, 0; y) 256 y^ H{0, 0, 1, 0; y) 



i^ + yy 



{i-y'f 



l+y\ il + y'fH{0;y) 



1 - r 



+ 



{i-y'Y 



l + 4y-4y3-^4_4^^_8^2^^_4^4^^ 

(i-y'f 

(_l_8y-8|/3 + ^4) ^(0;y) , 2{l+y^) H{l;y) 



{i-y' 



, 3(l+y2)2ff(o^O.^) 2_(1 + y2)- iJ(o, 1; ^) 



+ 



,2\2 



1-1/2 



(1-y^f 

4(l+y2)2i^(i^0;t/) 



(1-y 



2^2 



(1 - y^)^ 
1 



{l-yf{l+yf 



3-7y+10y' + 10y^-7y^-3y'-C2 

+41t/C2-18y'C2 + 18y'C2-45t/^C2 + 5?/'C2 + 8C3 
-24 y C3 - 56 Cs - 56 Cs - 24 Ca + 8 Cs 

4 (2+ Ij/- ljy^^-2jy^ + (. + 2jy2(, + ^i(,) H{-l:y) 

(1 - y^)' 

3 + 12 y + 36 y2 + 22 y^ - 47 y^ - 26 y^ 



(1-yf (i + yf L 

-10C2 + 6yC2-28y2C2 + 20y3C2-18y^C2 



+14y'C2 



^(0;y) 



(15 + 42y-42y3-15y^ + 16C2 + 32y2C2 + 16y^C2) H{l;y) 

(1 - y^)' 
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_2 (5 + ley + 6^ + 16^ + 57/^) H{-l,0;y) 

(1 + 31 1/ - 65 _ 9 y3 + 10 1/5) H{0, 0; 2/) 

^ (1-1/)' (1 + 2/)' 

2 (3 + 22y + 20y2 + 22y^ + 5|/^) g(0,l;|/) 

(5 - 58 y - 2 - 58 + 9 y^) //(I, 0; y) 

4 (1 + y') H(l, 1; y) 8 (1 + y^)' i/(-l, 0, -1; y) 

(1-y^)' 

8 (-1 + 3 y + 8 y^ + 12 y=^ + y^ + y5) H{-1, 0, 0; y) 

{i-yf{i + yf 

8 (1 + y2)' H(-l, 0, 1; y) 4(1 + y) //(O, -1, 0; y) 



{1-y^f 1-y 
4 (-l + 2y2 + 3y^) i/(0,0,-l;y) 

(1 - y'f 

(-7 + 7y - 16y2 + 52 27 + 27 y^) i/(0, 0, 0; y) 

(l-y)^^(l + y)^ 
_2 (l + 3y- 10?/ + io?/_ 7^/+ 11^5) ii"(0,0,l;y) 

(l-y)'(l + y)' 
2 (-3-5y-36y2-16/-15/ + 7y5) ff(0,l,0;y) 

(l-y)-^(l + y)^ 

4(l+/)'F(0,l,l;y) 12(l + /)'iJ(l,0,0;y) 

8(l+y2)'iJ(l,0,l;y) 16 (1 + /)' ff(l, 1, 0; y) _ 
(1-/)' (1-2/^)' 



11 (1+/) 

12 (1 - y2) ' 

(l + y2) (-67 + 67y2 + 18C2 + 18y^C2) 
36(l-y2)2 

^ (1 + /) F(-l;y) _ 2y2i/(0;y) _ (l + y^) i/(l;y) 
1 — y^ 1 _ y2 1 — y^ 

(1 + y^)' g(0, -1; y) 2 y^ (1 + /) (0, 0; y) 
(1-2/^)' (1-2/^)' 
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3 2 fl21 + 1133?/ - 1254?/^ - 1254 + 1133 y"" 

54 (1 - y) {l+y) L 

+ 121 - 540 1/2 (2 + 540 y^ C2 + 540 y^ C2 - 540 C2 
-351 Cs + 351 2/ Ca - 594 y^ (a - 702 y^ (3 + 405 y^ Ca 
-405^5(3 

2(i + y2) (i-y^ + (2 + ^^(2) 
(1 - y'f 

^37^^ + 74 y - 78 y^ - 24 + 11 - 50 

18(l-y) (l + y) L 

-18C2 + 18yC2-108y2C2 + 396y3C2 -234y4C2 

+234y5C2 
1 

^18(l-y2 

+72y2C2 + 36y^C2] H{l;y) 



H{0;y) 

- 53- 102y + 102y=^ + 53y^ + 36C2 



2 (1 + y2) //(-I, -1; y) 4 (-3 + y^) H{-1, 0; y) 



1 _ y2 1 — y^ 

6 (1 + y2) i/(-l, 1; y) 4y (1 - y + 2 y^) H{0, -1; y) 



1-?/' (l-y)'(l + y) 

(ll + 2y-53y2 + 74y=^ + 62y^) H{0,0;y) 
6{l-yf (l + y) 



3(1 


-yf (l + y) 


+ 1- 


- y^ 


(-61 + 175 y + 


35 y'^ + 79 y3) 


i/(l,0;y) 52(l + y2 


) i^(i,i;i/) 


6(1 


-y)' (l + y) 


3(1 




4{l + y'fH{- 


-1,0, -l;y) 


2 (-5-4y2 + y4) H{- 


1,0, 0;y) 



(1 - y^)' (1 ~ y'f 

^4^ + y'f H{-1, 0, 1; y) 2(1 + y'f H{0, -1, -1; y) 



(1-1/2)' (1-1/2)2 
2 (-3 + 3y + 8y2 + 16y^-y^ + y5) g(0,-l,0;y) 

(l-y)^(l + y)2 
6(l + y2)2g(0,-l,l;y) 2 (1 + 8y2 + 7y^) g(0, 0, -1; y) 

(1-1/2)2 (l-y2)2 
3y2 (-l + 7y-4y2 + 4y3) ff(0,0,0;y) 
(l-y)='(l + y)2 
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2 (-1 + y - 22 1/^ + 2 y-^ ~ 1 1 1/ + 11 1/5) (0, 0,1;^) 
6(l+|/^)^g(0,l,-l;2/) 

2 (-l+2/-2/2 + i5y3_7^4^7^5) i7(o,l,0;2/) 



(l-|/)-^(l + y)^ 



10(l + y2)2ij(o,l,l;j/) 4(l + y2)-i7(i,o,-l;|/) 



.2 I „,4 



+ 



{i~y^r 



2 (3 + 4|/2 + y4) iJ(l,0,0;2/) 4 (1 + y^)^ ^^(i^ q, 1; ?/) 



{i-y^y 



{i-y^r 



(82) 



1 + 



3(l-y2) ' 

5(1 + 1/^) 
9(l-?/2) ' 



4 (7 + 48|/ + 7y2) 2 (5 + 12y + 5y2) i7(0;7/) 



27 (1 - 



+ 



9 (l-y2) 



4 (5 + 12|/ + 5|/^) g(l;|/) 2(1 + y^) H{0, 0; y) 



9 (1 - y2) 



3 (1 - y') 



^ 4(l+y2) H{0,l;y) ^4(1 + y') H{1, 0; y) 



3 (1 - y2) 
8(l + y2) //(l,l;y) 



3 (1 - y') 



3 (1 - y') 



0; 
0; 



27(1-2/2) 



14 - 72 y + 106 + 384 + 106 y^ -72y^ + 14 y^ 



+216 2/2 C2 - 216 2/^ C2 + 81 y Cs - 324 y' - 378 y^ (s 
-3242/^C3 + 8l2/'C3 
, 16yH{-l;y) 



+ 



1-2/2 



9(i-#(i + y)^ 



5 - 232/ - 101 2/2 + 531 2/^ + 272/ - 3532/^ 



+229 y6 + 5 + 9 2/ C2 - 45 y2 C2 + 378 y^ (2 - 378 y^ (2 



+45 2/'C2-92/'C2 



i/(0;2/) 



32yH{l;y) 128 2/^ iJ(-l, 0; 2/) , 642/^ i7(0, -1; 2/) 



1-2/2 



{i-y'Y 

31 



+ 



(l-Z/2 



(83) 
(84) 



(85) 

(86) 
(87) 



2 (-1 + 481/2 + ^/4) H{0,Q;y) 128y' H{0, l;y) 



{1-yf (l + y) 
Sy {l-Ay-26y'-Ay^ + y^) H{0,-l,0;y) 

Sy {l-Qy + y^) g(0,0,-l;y) 

(l-yf (l + y) 
2 (-1 + y)y H{0, 0, 0; y) 256 y^ H{0, 0, 1; y) 



4.2.2 The Pseudoscalar Case 

KL2 = 3?6-2 ; (89) 

^ ^ (l + y^)g(0,0;y) ^ 2 (1 + y^) g(l, 0; y) 

1 _ y2 1 — 

3(l + y2)^i/(0,0,0;y) 2 (1 + y^)^ i/(0, 1, 0; y) 

(1 - y^)' (1 - y'f 

^4^1 + y^ f H{1. 0. 0: y) 12 (1 + r')' -^(1: ?/) C2 



(1 - y^)' (1 - y'f 

(1 + y^) //(O; y) (-1 + y^ + 10 C2 + 10 y^ C2) 

(1 - y'f 

-l + y^ + 4C2 + 4y^C2 
l-y2 

-2 (5 + 12y-18y^ + 12y=^ + 5y^) //(-l,0,0;y) 
4(l + y2) //(0,-l,0;y) 

(-1 - 9 y + 9 y^ + 71 y^ + 16 y^ + 10 y^) 1/(0, 0, 0; y) 

(l-y)^(l + y)=^ 
2 (3+14y-4y2 + 14y=^ + 5y^) ij"(0,l,0;y) 

ii-y^f 

(5-14y + 46y2-14y3 + 9y^) i/(l,0,0;y) 

(T^V? 
4 (l + y^) i/(l,l,0;y) 



l-y2 

8(l + y2)'if(-l,0,-l,0;y) 



(90) 



{i-y'f 
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8 {l-2y + 2y' + 2y^ + y^) 0, 0, 0; 

8(l + |/^)'ff(-l,0,l,0;y) 

4 {l + Ay + 2y^ + Ay^ + y^) g(0, -1, 0, 0; |/) 

{l-y) {l + yf 
4 {-l + 2y^ + ^y^) iJ(0, 0, -1, 0; 

(1-1/^)' 

(7 + 7y + 24y2 + 44y3 + 27y4 + 27 1/^) i/(0, 0, 0, 0; y) 



(i-yr(i + yr 

2 (-1 -5y + 10y2 + i0y3 + i5^4^;^i^5) 7/(0,0,1,0;^) 



(i-yr(i + yr 

2 (3 + lly + 12y2 + 8y=^-y4 + 7y5) //(0, 1, 0, 0; y) 



4(l + y2)'i/(0,l,l,0;y) 



(1 - y^)' 
12(l+y2)2//(i^0,0,0;y) 



(1 - y^)' 

8(l + y2)'i/(l,0,l,0;y) 



(1 - y^f 

16(l + y2)'//(1.1.0.0:y) 



2 (-5 + 25y-17y^ + 13y=^) //(l;y)C2 

(l-y) (l + yf 
12 (l + 4y + 2y^ + 4y3 + y4) g(0,-l;y)C2 

(l-y)(l + y)' 
8 (4 + 10y + lly^ + 5y=^-2y4 + 4y5) ff(0,l;y)C2 

(i-y)'(i + y)' 

48 (l + y2)2i/(i, 1.^)^2 



(i-y'f 



^4i7(-l,0;y) 



2-2y^ + 5C2-8yC2 + 10y'C2 



+8^3(2 + 5y^C2 
//(l,0;y) (15 -2y + 2y3_ 15^4 ^40^^ ^80^2^^ ^40^4^^) 



i/(0,0;y) 



(l-y)^(l + y)^ 



(1-y^f 

-3-14y + 6y2 + 8y3-35y^-26y^ + 24C2 
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+28 y C2 + 76 (2 + 108 (2 + 64 1/^ (2 + 68 1/^ C2 



2H{-l-y) 



(1 - y^) 

-2C3-4y2C3-2y^C3 
1 



15 C2 + 36 1/ C2 - 54 C2 + 36 C2 + 15 1/^ C2 



+ 



145 + 145 7/ - 290 7/^ - 290 + 145 + 145 



20(l-yr(l + yr 
-160 C2 - 600 y C2 + I960 C2 + 2120 (2 - 1800 (2 
-1520 y^C2 - 480 y ln(2) C2 - 1440 ln(2) C2 - 1440 ln(2) C2 
-480 y^ ln(2) C2 + 272 C2' + 376 y C2' + 912 y^ (2' + 416 y=^ C2' 
+288 y^ C2' + 392 y^ C2' - 280 Cs - 280 y' Ca - 280 y^ Cs 
+80y^C3 - 200y5C3" 



H{0;y) 



(i-yr(i + yr 



3 + y-2y2 + 2y3-y4-3y5-C2-21yC2 



+16y2C2 + 140y3C2 + 33y"C2 + 25y^C2-8C3 + 8yC3 

+8y\3-8y\3-8y^C3 + 8y^C3 



3?c_2 ; 



3?c_i ; 

-2 (l + y^) if(-l,-l,0;y) 
l-y2 

4 (3 + 2y-4y2 + 2y3 + y^) ff(-l,0,0;y) 

6 (l + y^) g(-l,l,0;y) 
l-y2 

4y(l-y + 2y2) g(0,-l,0;y) 

(1-yf (1 + y) 
(-11 - 35 y - 45 y^ + 189 y^ + 128 y^ + 62 y^) //(O, 0, 0; y) 
^ 6(l-yf(l + yf 

2 (-l + 15y-15y2 + 15y=^ + 34y^) H{0,l,0;y) 

3 (1 - y^)' 
6 (l + y2) i7(l,-l,0;y) 

l-y2 

(-61 + 6 y + 66 y2 + 6 + 79 ^(^^^ q, 0; y) 

^ 6(l-y2f 
52 (l+y2) g(l,l,0;y) 
3(l-y2) 
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4(l + y^)^g(-l,0,-l,0;y) 

(1 - y'^f 

2 (-5-4|/^ + j/^) if(-l,0,0,0;y) 

(1-1/^)' 
4(l + |/^)^g(-l,0,l,0;y) 

2(l + y^)^g(0,-l,-l,0;y) 
(1-1/^)' 

2 (3 + 3|/ + 81/^ + 1/^ + 1/5) i/(0,-l,0,0;|/) 

(1 + l/f 
6(l + y^)'iJ(0,-l,l,0;i/) 

2 (l + 8t/^ + 7|/^) g(0,0,-l,0;y) 
(1-1/^)' 

1/^ (7+17|/+12|/^ + 12|/3) iJ(0,0,0,0;|/) 

2 (l + y+14y^ + 10y=^ + lly^ + lly5) //(O, 0, 1, 0; y) 

6 (1 + y^f g(0, 1,-1,0; y) 

2 (1 + + 5 j/^ + 11 j/^^ + 7 j/^ + 7 //((). 1.0. 0: i)) 
10(l + j/^)'f/(0.1.1.0:y) 

(1 - y^f 

4(l + jy^)-//(1.0.-L0;y) 

2 (3 + 4y^ + y^) //(I, 0, 0, 0; y) 

4(l + y^)^//(l,0,l,0;y) 

(1 - y"? 

(19 + 19y + 2y^ + 50y3 + 7y" + 7y^) //(0,-l;y)C2 

(l-yf (l + yf 
(11 + 11 y + 34 y^ + 58 y^ + 35 + 35 y^) ij-(o, i; y) 

(1-yf (1+l/f 
2 (1 + y^) i/(-l,0;?/) (-l+y^-llC2 + Z/^C2) 
(l-l/^)' 
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+ 



H{l,0;y) 

^ H{0,0;y) 
18{l-yf{l + yf 

+18 y C2 + 504 C2 + 864 y^ C2 + 666 y^ C2 + 666 C2 

2 (l + y2) H{l;y) (-52(2 - 9y C2 + Qly\2 + ^Cs + ^y'' Cs) 



- 53 + 54y-54y3 + 53y^ + 252(2 + 360^2(2 



-67- 148y-36y2-90y^- 185y^-50y^ + 18C2 



+ 



H{-l;y) 

-2C3-4y2C3-2y^C3 
1 



3(i-yT 

37 C2 + 24 y C2 - 48 y2 C2 + 24 y3 C2 + 1 1 2/' C2 



- 4345 - 4345 y + 8690 y^ + 8690 y^ - 4345 y^ 



540(l-y)^(l + yy 
-4345 y^ + 8760 C2 + 12000 y C2 - 21960 y^ C2 - 15480 y^ C2 
+13200 y^ C2 + 3480 y^ C2 + 6480 y ln(2) C2 + 19440 y' ln(2) C2 
+19440 y^ lii(2) C2 + 6480 y^ lii(2) C2 - 1701 C2' - 1701 y C2' 
-16902 y2C2^ - 13014 y^C2^ - 13257 y^ (2^ - 13257 y^ (2^ 
+4410 C3 + 4950 y Ca + 2700 y' Cs + 2700 y' (3 - 630 y^ (3 
-1170 y^Ca 

H{0;y) 
54(l-y)2(l + yr 



121 + 283 y + 162 y^ - 162 y^ - 283 y^ - 121 y^ 



-198 C2 - 630 y C2 - 1350 y' C2 + 2862 y^ (2 + 2844 y^ C2 
+1656 y' C2 - 351 Ca - 351 y Ca - 162 y^ Ca + 270 y^ Ca 

+405y^C3 + 405y^C3 



(94) 



3?(i_2 



10 (1 + y2) H{0, Oyy) ^ 20 (1 + y^) //(I, 0; y) 



9(l-y^ 



9 (1 - y^ 



2 (1 + y^) HjOMy) , 1(1 + y') H{0, 1, 0; y) 



3(l-y2 



3 (l-y2 



^ 4 (1 + y^) H{l,0,0;y) ^ 8_(1 + y^) i7(l, 1, 0; y) 



3(l-y2) 3(l-y2) 
16 (1 + y2) H{1; y) C2 4(1 + y^) i7(0; y) (-7 + 9 C2) 



3(l-y2) 



27 (l-y2) 



(95) 
(96) 
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-49 + 49 + 84 C2 + 156 C2 + 36 Ca + 36 y^ (s 



27(1-1/2) 



3?e_2 = ; 



3?e_i = ; 



2 (l + 2y + 26y2 + 2y3 + ^4) H{ 0,OAy) , By //(-I, 0, 0, 0; y) 



3(i-y) (1 + y)' ' 

8y/7(0,-l,0,0;y) 8y i/(0, 0, -1, 0; y) 12 y i7(0, 0, 0, -1; y) 



1 — 1/2 1 — y^ 1 ~ 

2yi/(0,0,0,0;y) 16y i/(-l, 0; y) (2 24y if(0, -1; y) C2 



l-y2 
2i/(0,0;y) 



9(l-y)^(l + y)^ 



\ — 1 — y'^ 

5-22y + 59y2- 116y^ + 59y^-22y^ + 5y^ 



-9yC2 + 18y'C2-18y'C2 + 9y^C2 
4i/(0;y) 



27(l-y2 



- 14 + 12y + 38y^ + 24y'' + 38y^ + 12y^- My*' 



+9C2 + 207y2C2 - 432y3C2 + 207y^C2 + 9y^C2-81yC3 

+162y^C3-81y^C3 
1 



+ 



135 (1 - y2) 



2035 - 120 y - 4310 y^ - 120 y^ + 2035 y* - 1260 C2 



+540 y C2 + 1440 y' (2 + 540 y^ (2 - 1260 y^ (2 - 783 y (2' 
+783 y^ (2'^ 



(1+y^) if(0;y) 2 (1 + y2) (1; y) 3 (1 + y^)^ iJ(0, 0; y) 
l-y2 + l-y2 + (i_^2)2 

2(l + y2)^if(0,l;y) 4 (1 + y2)2 ^ (1^ p. ^) 

(l-y2)^ (l-y2)2 
(1 + y^) (-l + y2 + 4C2 + 4y^C2) 



-2 (5+ 12y- 18y2 + 12y3 + 5y^) H{-l,0;y) 

{i-y'f 

4 (l + y2) iJ(0,-l;y) 
(1-1/)^ 
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_{-l-9y + 9y^ + 71y^ + 16 ?/ + 10 y^) H{0, 0; y) 
{1-yfil + yf 

2 {3+Uy-Ay^ + Uy^ + 5y^) H{0,l;y) 
(5- 14 7/ + 46 7/2- 14 7/3 + 9 1/4) H{l,0;y) 

4(1 + 1/2) H{l,l;y) 
l-y2 

8(l + y2)2ff(_i^0^_l.^) 

8 (l-2y + 2y2 + 2y3 + ^4) ff(_i^o^o.^) 

(1 - y'f 
8{i + y'fH(-i,o,-i-;y) 
(1 - y'f 

4 (l + 4y + 2y2 + 4y3 + ^4) ff(o^_i^o.^) 

(l-y) {1 + yf 
4 (-1 + 2 7/2 + 3/) g(0,0,-l;y) 

(i-y'f 

(7 + 7 7/ + 24 7/2 + 44 + 27 / + 27 y^) H{0, 0, 0; y) 

^ (l-?/)^(l+7/)-^ 

2 (-l-5y+ 107/2 + 107/''^ + 157/^ + 117/5) i/(0,0,l;7/) 

^ (l-y)'(l + y)' 

2 (3+ll7/ + 127/2 + 87/3-7/^ + 77/5) g(0,l,0;7/) 

(i-yf + 

4 (1 + 7/2)^^(0,1,1; 7/) 

12(1+7/2)2^(1^0^0-^) 

{i-y'f 

8(l + 7/2)^g(l,0,l;7/) 

(1-1/^)' 

16(l + 7/2)2ff(i^l^0;7/) 

(1 - y^)' 

4(1 + 7/2) g(-l:y) (2-27/2 + C2 + |/^C2) 

{i-y'f 

-7?^ ( - 3 - 17y^ - 9t/^ - 26/ + IOC2 + 47/C2 

(1-7/2) V 
+247/2(2 -4/C2 + I4/C2 
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H{l;y) (15-2y + 2|/=^-15y^+16C2 + 32|/^C2 + 16|/^C2) 

(1 - y'f 



(l-y) {1 + yf 

+4y^C2 + 6y'C2 + 5y\2 + 8Cs-8y\s + 8y^Cs\ ; (103) 
9c_2 = 9c_2 ; (104) 

^c_i = »c_i ; (105) 
-2 (1 + 1/2) i7(-l,-l;|/) 



3^ Co 



4 (3 + 2y-4y2 + 2y3 + ^4) ff(-i,0;^) 

6(l + y2) //(-l,l;y) 

Ay{l-y + 2y') g(0,-l;y) 

(l-y)' (1 + y) 
(-11 - 35 y - 45 + 189 ?/ + 128 + 62 y^) H{0, 0; y) 

^ 6(l-y)'(l + y)' 

^ 2 (-1 + 15|/- 15 y' + 15 1/3 + 34 1/4) i/(0, 1; y) 



3 (1 - y2)2 
6 jl + y^) H{l,-l;y) 

1 - 

(-61 + 6 y + 66 1/2 + 6 y3 + 79 y'l) if (i^ Q; |/) 

6(1-1/2)2 
52 (1 + 1/2) H{l,l-y) 

3 (1 - y') 

4(l + y2)2ff(_1^0^_l.^) 

(1 - y'f 

2 (-5- ly2 + yi) //(_LO^O;i/) 



(i-y'f 



4(l + y2)2ff(_1^0^1.^) 

(i-y'f 

2(l + y2)2ff(o^_l^_l.^) 

(1 - y'f 

2 (3 + 3y+8y3 + y^ + y^) H{0,-l,0;y) 

a-yfa + yf 
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6(l + |/^)^g(0,-l,l;y) 

2 (l + 8|/^ + 7|/^) g(0,0,-l;y) 

(7+17|/ + 12|/=^ + 12y3) if(0,0,0;i/) 
2 (1 + 1/ + 14 1/2 + 10 + 11^4^11 1/5) if(o,0, 1;|/) 



(1-^^(1 + 2/^ 



6(l + i/2)^g(0,l,-l;|/) 
2 (l + y + 5y2 + iiy3 + 7^4^7^5) i7(o,l,0;|/) 



(i-?/r(i+yr 



10 (1 + 1/2)^^(0,1,1; I/) 

(1 - Z/^)' 
4(1 + 1/2)2 i/(i^O^_l.^) 

(1-2/^)' 
2 (3 + 4y2 + y4) (1,0,0;^) 

(1 - V"? 
4(l + y2)2ff(i^0^1.^) 

(1 - y^f 

2 (l + y2) //(-l:y) (l-y2 + (, + y2^) 



18(l-y2)^ 
+36y^C2" 

, ^^(0;i/) 



18(1-7/2) 

+252y2(^2 + 234|/^C2 
1 



(1 - y^f 

53 + 54 1/ - 54 1/3 + 53 1/^ + 36 C2 + 72 C2 



- 67- 81y + 45 1/2- 135^3 _ 50^4^ 18^2 



121 - 283 y - 162 + 162 y^ + 283 



54(l-|/)'(l + l/)3 
+121 y^ + 540 1/2 (2 + 540 1/3 (2 - 540 y^ C2 - 540 y^ C2 
+351 C3 + 351 y Cs + 162 1/2 (3 - 270 y^ (3 _ 405 y^ Cs 

-4051/5(31 ; 



(106) 



; 



(107) 



(108) 
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10 jl + y') H{0;y) 20 (1 + y') H{1; y) 



28 (1 + ?/^ 
27(1-1/2 

2 (1 + y^) g(0,0;i/) 4 (1 



9 1 



y 



y 



9 (1-1/2) 
2^ //(0,l;i/) 



3 1 



r 



4 (l + y2) if(l,0;i/) 



3 (1 - 2/2 



3(1 
1+y 



-yv 

2^ //(1,1;2/) 



3(1-2/2 



(109) 
(110) 



53 e_i 



3 e_i 



:iiii 



3 Co 



2 (1 + 22/ + 26?/2 + 27/3 + 2/4) g (0,0;2/) ^ 8 2/g(-l, 0, 0; 2/) 



3 (1-2/) (1+2/)' 
8 2/i/(0,-l,0;2/) 8 2/ //(O, 0,-1; 2/) 



1-2/2 
2H{0;y) 
9{l-y)\l+yf 

+92/C2-182/'C 
4 

^27(1-2/)' (1 + 2/) 
-162 2/'C3 + 8l2/'C3 



1 - 
22 2/ 



r 



+ 



1-2/2 

2 2/^(0,0,0; 2/) 
1-2/2 



59 y^ - 116 2/^ + 59 2/ - 22 2/^ + 5 2/^ 



182/'C2-92/'C2 



14 - 40 2/ + 28 2/ - 40 2/^ + 14 2/^ + 812/(3 



;il2) 



The light quark loop insertions in Fig.|2(h) were studied before by [18,19], using a small quark 
mass for regularizing the double- and single logarithmic mass divergences of these contributions. 
As we use dimensional regularization setting the light quark masses to zero, only the coefficient 
of the leading singularity of the CfTrNj term in the scalar and pseudoscalar form factors, in 
Eqs. (71-73), (83-85) and Eqs. (95-97), (107-109) respectively, can be compared, and we agree 
with [18,19]. 

The scalar and pseudoscalar CpTji coefficients, respectively in Eqs. (74-76), (86-88) and 
Eqs. (98-100), (110-112), are obtained by adding the contributions coming from the diagrams 
in Fig. 121 (f) and (g). Our results for the latter, with the heavy quark loop insertion, agree 
with the ones in [19], if one takes into account that the results in [19] are listed for an on-shell 
renormalization of the coupling constant, while we use MS-renormalization for the coupling 
constant. 



5 Threshold Expansions 

In this section we provide the expansions of our results in the threshold limit s ~ 4m2 (2/ — > 1 
in the transformed variable). We define 




(113) 
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as the small expansion parameter and give the results up to the zeroth order 



5.1 One- loop Scalar Case 



3?a_i = ; 



3?ao 
3?ai = 2 



-1 + 



3C2 



/3 ' 

6C2 (-l + ln(2)+ln(/3)) 
/3 



a_i 



1 

~2^ ' 

-1 + ln(2) + ln(/?) 



1 



/3 



-2 + 2 ln(2) - ln'(2) + C2 + 2 ln(/3) 



-2 ln(2) ln(/5) -ln2(/3) 



5.2 One- loop Pseudoscalar Case 



3fJa_i = 3?a_i ; 



3fJ oo 



-3 + 



3C2 



2- 



/3 ' 

6(ln(2)+ln(/3))C2 

P 



a_i 



9ao 



$5 Oi 



ln(2) + ln(/3) 

-ln^(2) - 2 ln(2) ln(/3) - ln^(/3) + C2 
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Two-loop Scalar Case 



3? 60 



3C2 3C2 



3(2 (-l + ln(2) + ln(/3)) , 3 (2 (1 + 4 ln(2) + 4 ln(/3)) 



3C2 



/32 + 2 

-4 + 8 ln(2) - 41n2(2) + 3(2 + 8 ln(/3) 



2/32 

-8 ln(2) ln(/?) -41n2(/?) 

3(2 



/3 



+ 



5 + 98 (2 - 40 ln(2) (2 - 48 W(2) (2 + 36 (2' 
-44 (3 - 88 (2 ln(/3) - 96 ln(2) (2 ln(/3) 
-48(2ln2(/3)" 



KC_2 

3? Co 



0; 



0; 



(2 (-97 + 66 1n(2) + 66 ln(/3)) 
6/3 

+|-tl^±4M?))^-5C3-4C..u(«; 



3f?(i_2 = ; 



3f?ci_i = ; 



2(2 (-ll + 6 1n(2) + 6 1n(/3)) 5 



3/3 



9 ' 



3?e_2 = ; 



3?e_i = ; 



5Reo = li^-^ + 41n(2) (-1 + 4(2) ; 
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2/52 2p 

6C2 (-l + ln(2)+ln(/3)) ln(2) + ln(/3) 
+C2 (11 + 12 ln(2) + 12 ln(/3)) ; 



$5 Co 



11 

24^ ' 
_ 

~72^ 
1 



- 239 + 291 ln(2) - 991n^(2) + 291 ln(/3) 



54 /?L 

-198 ln(2) ln(/3) -99W(/3) 

+2C2 ; 



1 

5 



18/3 ' 
2 



31 - 33 ln(2) + 9111^(2) - 33 ln(/3) 



27/3 

+18 1n(2) ln(/3)+91n2(/3)| ; 



$>e_2 = 0; 



$Je_i = 0; 



$560 = 2 



Two-loop Pseudoscalar Case 



3?6_2 = 3?6-2 ; 



3 (ln(2) + ln(/3)) (2 3 (1 + 4 ln(2) + 4 ln(/3)) C2 



/32 



+ 



44 



-9(2 , 3C2 (-41n^(2)-8 1n(2)ln(/3)-41n^(/3) + 3C2) 

= -^—^ 77^ 



1 

+ 4 



29 - 26 C2 - 24 ln(2) (2 - 48 ln\2) (2 - 120 ln(/3) (2 
96 lii(2) ln(/3) C2 - 48 In^P) C2 + 36 C2' - 72 C3] ; (152) 



3?c_2 = 0; (153) 

3?c_i = 0; (154) 

m = -g-izlI±Z|M^_i2 1n(,)C2 

(-31 + 66 1n(2) + 66Jn(^_^2C3; (155) 



6(3 

3? J_2 = ; (156) 

3?J_i = 0; (157) 

a;j„ ^ 1 2(-5 + 6 1n(2) + 6 1n(/^))C.^ (^^3, 

3fJe_2 = ; (159) 

3fie_i = ; (160) 

^eo = y-3C2 + 12 1n(2)C2-^ ; (161) 

$5 6_2 = ; (162) 

^l-. = ^-3C.-||; (163) 

_ -1 - 3 ln(2) - 3 ln(/3) , 6 (ln(2) + ln(/3)) (2 
^^0 - ^ + ^5 

+3 (5 + 4 ln(2) + 4 ln(/3)) C2 ; (164) 

= Q'c_2 ; (165) 
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(166) 



3^ Co 



-47 + 93 ln(2) - 99 ln^{2) + 93 ln(/?) 



54/3 

-198 ln(2) ln(/3) -991n^(/3) 



+ 6C2 ; 



Qdo -- 



(167) 
(168) 
(169) 



J-15 1n(2) + 91n^(2)-15 ln(/3) 
+18 ln(2) ln(/3) + 91n2(/3) 



»e_2 = ; 
Qe-i = ; 
^^eo = 3C2 • 



(170) 
(171) 
(172) 
(173) 



6 Asymptotic Expansions 



In this section we provide the expansions of the form factors in the hmit s ^ (t/ ^ in the 
transformed variable). We define r — s/m^ and L = In (r) and keep terms up to the second 
order in 1/r. 



6.1 One- loop Scalar Case 

3fJ ao = 



-3 + 2L 2 ^ ^ 
1 + L ; 



19 + 32 L - 2 + 16 C2 



;i74) 



2r2 

^ 2 (1 + 4L) ^ -2-L^ + 8C2 



3? ai = 



r 2 
3 + 456L-96L^ + 4L3 + 768C2-96LC2 + 48C3 

r 6 

+2C3; 



(175) 



(176) 
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2 



$5 01 



-38 + 161-12 + 4(2 

^ 2 (-3 + 4L) ^ -2-L^ + 4C2 



One-loop Pseudoscalar Case 



Sa. = 2(l + 2L)^-2-L- + 8C. 



r 2 



-3 + 16L - 2L2 + 16C; 



2r2 

^-a, = -2 + L + ^-^i^^±^^^^-4LC2 + 2C3 
r 

99 + 72L-48L^ + 4L^ + 384C2 -96LC2 + 48C3 

^ 12^2 



Q^Oo = ; 
'^tti = ; 



Two-loop Scalar Case 

5-5L + 2L2-12C2 -2 (-1 + L) 

0-2 = t; 1 



-2 



1 - 2 L + - 6 C2 
+ 2 ' 

17-87L + 37L2-4L3- 232C2 + 80LC2 



-6 L + 9 L2 - 56 C: 



2r2 
2 



+ 

r 
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20235 + 6240 L + 10080 - 1760 + 150 

120 L 

-58680 C2 + 41760 L C2 - 7200 C2 - 34560 ln(2) C2 
+10224 C2^ + 10320 C3 - 10560 L Cs 



1 



30 - 360 L - 125 + 480 C2 + 2580 L C2 - 30 C2 



15 r 

-1080 ln(2) C2 + 78 (2^ + 210 Cs - 240 L C3 
29 L3 7L4 „.^68C2' 3L^ (-I + 8C2) 

-uCs + L (-3 + C2 + 8C3) ; 



11 (-3 + 2 L) ^ 11 -11 (-1 + L) 



12 r2 6r 12 

-165 + 188 L - 36 + 12 + 234 (2 - 216 L C2 - 72 (3 



+ 



36 

-67+I8C2 



18 r 

-49 + 67L+18C2- I8LC2- I8C3 
36 

1 



58815 + 264770 L - 20010 + 2580 - 270 



2160 r2 L 

-23520 C2 - 63360LC2 + 151201.2(2 + 311040 ln(2) C2 
-55296 C2^ + 155160 (3 - 51840 L (3 



1 



628 + 2274 L - 261 + 126 C2 - 108 L C2 



54 r 

+ 1944 ln(2)C2 + 1296 Cs 

+4- \ - 4345 + 1210 L - 1005 + 165 + 8760 C2 
540 L 

-1980 L C2 + 270 C2 - 1701 (2^ + 4410 (3 
-3510 LCs 



-3 + 2L 2 ^ -1 + L 



3r^ 3r 3 
5 (-3 + 2 L) ^ 10 ^ -5 (-1 + L) 



Qj.2 9^ 9 



48 



477 - 2224 L + 348 - 12 - 2568 C2 + 144 L (2 - 144 (3 

54 

8 (-5-39L + 9L2-63C2) 
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7 Summary and Outlook 



If Higgs boson will be discovered, then the next important task will be the experimental determi- 
nation of its couplings to the other known fundamental particles, in particular to heavy quarks 
Q = c,b, t. Therefore precise predictions of Higgs boson decays into heavy quark- ant iquark 
pairs are mandatory, not only for the partial decay rates, but also for differential distributions. 

To put forward this aim we have computed the two-loop QCD corrections to the decay 
amplitude of a neutral Higgs boson, with scalar and pseudoscalar couplings to fermions, into 
a heavy quark- ant iquark pair, h — > QQ. This amplitude is determined by two form factors, 
which we calculated with the algebraic and analytic techniques employed in [35,36,39,44,45], for 
arbitrary quark mass m and squared four- momentum q'^ of the spin-zero particle. These scalar 
and pseudoscalar vertex functions have also a number of other applications, including resonant 
Higgs boson production in high energetic muon collisions, ii'^n" — > QQX. Applications to 
hadron physics include matrix elements of scalar and pseudoscalar quark currents, serving as 
interpolating fields for heavy mesons, between heavy (anti)quark states. Moreover, from the 
results presented here one obtains as a special case also the second order QED corrections for 
the amplitudes h — > //, / being a quark or a lepton. 

In addition, we expanded the scalar and pseudoscalar form factors near the QQ threshold 
and in the regime in^/q^ <^ 1. The former expansions may be applied, for instance, to the 
decay of heavy nonstandard Higgs bosons into tt, the latter ones apply to, e.g., h bb, cc. 

The renormalized form factors were expressed in terms of the MS coupling as for {Nf + 
1) quark flavors and the on-shell mass m of the heavy quark. They may be expressed in 
straightforward fashion as functions of the running MS mass m(A*), too. It is expected that 
the use of this mass parameter leads to a faster decrease of the coefficients of the perturbation 
series [4]. However, this is no general rule (see, e.g., [14]); the size of the coefficients of the 
series when using m or m must be studied for each observable separately. 

The renormalized form factors contain soft and coUinear divergences, which in our IR- 
regularization scheme appear as poles in e. These divergences have to be canceled against the 
divergences arising from the real radiation of gluons and massless quark-antiquark pairs, which 
in this paper were not taken into account. The IR singularity structure of QCD amplitudes 
involving massive quarks is so far not understood beyond the 1-loop level. The presently 
known two-loop heavy-quark vertex functions may provide insight into the generic singularity 
structure of the two-loop QCD amplitudes with massive quarks. As to the two-loop scalar 
and pseudoscalar form factors given above we notice that the infrared-divergent contributions 
proportional to the colour factor Cp show Yenni-Frautschi-Suura exponentiation [50], which 
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is characteristic of Abelian gauge theories. This is completely analogous to the case of the 
two- loop QCD vector and axial vector vertex functions [36]. However the IR-divergent terms 
proportional to the other colour factors cannot be explained by naive Abelian exponentiation. 
It remains an open problem whether an infrared factorization formula (as derived for massless 
QCD in [51,52]) exists for the generation of these singularities beyond the one-loop level. 

As emphasized, the results of this paper are one building block towards a differential de- 
scription of Higgs boson decays, h —> QQX, to second order QCD. 
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